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The. Rotorcraft  Flight  S.imubtion  is  a  multidisciplinary  mathematical  model  that  may  be  used  to  simu¬ 
late  a  wide  variety  of  helicopter  or  V/SfOL  aircraft  configurations  using  a  digital  computer.  Aircraft 
performance,  stability  and  control,  and  maneuver  characteristics,  as  well  as  rotor  blade  loads,  may  be 
estimated  using  this  analysis.  The  fuselage,  main  rotor,  tail  rotor,  wing,  elevator,  fin/rudder,  jet  thrust, 
and  weapon  recoil  are  treated  as  separate  aircraft  components,  allowing  detailed  representation  of  the 
aircraft  for  design  or  detailed  analysis  applications.  Six  rigid-body  fuselage  degrees  of  freedom  and  up 
to  six  rotor  blade  elastic  degrees  of  freedom  for  each  of  two  rotors  arc  accounted  for. 

Input  for  the  simulation  is  divided  into  logical  blocks  in  an  casy-to-undcrstand  format.  The  rotor  blade 
elastic  degrees  of  freedom  arc  omitted  if  stiffness  and  mass  properties  are  not  known.  Output  includes 
aircraft  trim  attitude,  control  positions,  performance  rotor  loads,  stability  and  control  characteristics, 
and  detailed  maneuver  response. 

While  correlation  of  the  predicted  results  from  the  Rotorcraft  Flight  Simulation  is  generally  satisfactory, 
especially  for  single-rotor  helicopters,  no  detailed  comparison  with  blade  loads  test  data  has  been  made. 
Current  development  of  the  analysis  includes  such  a  correlation  as  well  as  an  investigation  of  advanced 
control  systems,  development  of  an  improved  unsteady  aerodynamic  representation,  the  inclusion  of  a 
more  detailed  aerodynamic  representation  of  all  aircraft  components,  and  the  implementation  of  ad¬ 
vanced  numerical  techniques  that  will  result  in  reduced  run  times. 
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SUMMARY 

The  Rotorcraft  Flight  Simulation  Computer  Program  C-81  is  a  multidisciplinary  mathemati¬ 
cal  model  that  may  be  used  to  simulate  a  wide  variety  of  helicopter  or  V/STOL  aircraft 
configurations  using  a  digital  computer-  Aircraft  performance,  stability  and  control,  and 
maneuver  characteristics,  as  well  as  rotor  blade  loads,  may  be  estimated  using  this  model. 
The  fuselage,  main  rotor,  tail  rotor,  wing,  elevator,  fin/rudder,  jet  thrust,  and  weapon 
recoi!  are  treated  as  separate  aircraft  components,  allowing  detailed  representation  of  the 
aircraft  for  design  or  detailed  analysis  applications.  Six  rigid-body  fuselage  degrees  of 
freedom  and  up  to  six  rotor  blade  elastic  degrees  of  freedom  for  each  of  two  rotors  are 
accounted  for. 


Input  for  the  simulation  is  divided  into  logical  blocks  in  an  casy-to-undcrstand  format. 
The  rotor  blade  clastic  degrees  of  freedom  arc  omitted  if  stiffness  and  mass  properties 
are  not  known.  Output  includes  aircraft  trim  attitude,  control  positions,  performance, 
rotor  loads,  stability  and  control  characteristics,  and  detailed  maneuver  response. 

While  correlaf'jn  of  the  predicted  results  from  the  Rotorcraft  Flight  Simulation  is 
generally  satisfactory,  especially  for  single-rotor  helicopters,  no  detailed  comparison  with 
blade  loads  test  data  has  been  made.  Current  development  of  the  analysis  includes  such 
a  correlation  as  well  as  an  investigation  of  advanced  control  systems,  development  of  an 
improved  unsteady  aerodynamic  representation,  the  inclusion  of  a  more  detailed  aero¬ 
dynamic  representation  of  all  aircraft  components,  and  the  implementation  of  advanced 
numerical  techniques  that  will  result  in  reduced  run  tim^. 
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INTRODUCTION 


The  Rotorcraft  Flight  Simulation  Computer  Program  C-81  is  a  highly  versatile  digital 
flight  simulator  that  may  be  used  to  study  steady-state  and  maneuver  performance,  rotor 
loads,  stability  and  control,  aircraft  attitude  and  flight  path,  gust  response,  and  many 
other  aircraft  characteristics.  This  report  describes,  without  exhaustive  analytical  detail, 
the  history,  content,  use,  and  current  development  of  the  Rotorcraft  Flight  Simulation 
program.  The  analysis  and  use  of  the  program  are  further  explained  in  Reference.1;  1 
through  7. 
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HISTORY 


The  C-81  family  of  helicopter  flight  simulation  programs  has  been  under  development  by 
Bell  Helicopter  Company  and  the  Government  for  the  past  decade  (Figure  1).  Through¬ 
out  this  development,  certain  guidelines  have  been  followed.  First,  the  analysis  is  suffi¬ 
ciently  general  to  describe  a  wide  variety  of  helicopter  configurations  (conventional, 
compound,  tandem,  side-by-side)  for  a  bread  range  of  flight  conditions  (hover,  transition, 
cruise,  high  speed,  maneuver).  The  overall  analysis  has  a  uniform  texture;  i.c.,  the  level 
of  complexity  of  the  different  phases  (aerodynamic,  dynamic,  rotor  analysis,  etc.)  is 
uniform.  The  program  is  applicable  to  diverse  types  of  analysis  such  as  performance, 
stability  and  control,  and  rotor  loads.  Finally,  the  program  is  user-oriented  in  terms  of 
preparing  the  input  data  and  interpreting  the  results. 

The  first  major  step  in  computerized  analysis  was  a  digital  pro-am  to  determine  the 
overall  helicopter  performance  and  rotor  blade  bending  moments  for  level-flight  conditions. 
Aerodynamic  considerations  included  compressibility,  stall,  and  reversed-flow  effects  from 
two-dimensional  airfoil  test;  separate  treatment  of  specified  radial  segments  with  special 
lift  and  drag  characteristics;  and  procedures  to  balance  all  forces  and  moments  to  satisfy 
the  requirements  of  trimmed  flight.  The  introduction  of  coupling  between  in-plane  and 
out-of-planc  blade  deflections  in  the  roto.  dynamic  analysis  significantly  improved  the 
calculation  of  natural  frequencies  and  force '  response  for  rotor  systems  with  several  com¬ 
binations  of  number  of  bit.  les  *nd  types  of  hub  construction. 

The  next  major  development  of  the  analytical  simulation  included  the  addition  of  a 
maneuver  capability  with  six  degrees  of  freedom  for  the  helicopter  fuselage.  Definition 
of  the  airframe  was  extended  to  include  physical  dimensions  for  eg  location,  mast  length, 
and  tilt;  and  the  sizes  and  locations  of  wings,  elevator,  vertical  fin,  and  pylon  fairing. 
Contributions  to  lift,  drag,  and  side  forces  and  to  pitching,  yawing,  and  rolling  moments 
were  treated  separately  for  each  aerodynamic  surface  to  obtain  a  useful  method  of  cal¬ 
culating  stability  derivatives  and  maneuver  capability.  Control  linkage  ratios,  engine  power 
controls,  and  external  gust  disturbances  were  added  to  simulate  a  wide  variety  of  VTOL 
maneuvers. 

Under  the  contract  reported  in  Reference  2,  the  simulation  was  further  expanded  to 
encompass  ail  of  the  basic  rotorcraft  configurations:  single  main  rotor  plus  antitorque 
tail  rotor,  tandem,  sidc-by-sidc  tilting  rotor,  and  coaxial.  The  detailed  aerodynamic  and 
dynamic  treatment  of  the  second  rotor,  plus  provisions  for  locating,  orienting,  and  con¬ 
trolling  both  of  the  rotors,  led  to  an  all-purpcsc,  generalized  analytical  ability.  Two-, 
three-,  and  four-bladed  rotors  were  considered  for  hub  types  that  were  either  teetering, 
gimbaled,  articulated,  or  rigid  (hingelcss).  The  effects  of  gradual  penetration  of  a  shaped 
gust  field  by  the  rotor  disc  were  added  to  the  analysis  and  evaluated  during  the  course  of 
the  study. 

»n  1969,  under  an  Air  Force  contract,  a  version  of  the  analysis  was  prepared  with  a 
special  treatment  of  the  rotor  dynamics  to  allow  the  study  of  slowed-  and  stopped-rotor 
VTOL  configurations.  A  time-variant  analysis  was  added  that  accounts  for  the  rigid-blade 


flapping  motions  of  up  to,  and  including,  seven  blades  while  simulating  teetering, 
gimbaled,  articulated,  or  rigid  hub  configurations. 

A  time-variant  aeroclastic  rotor  response  analysis  based  on  the  modal  technique  was  in¬ 
corporated  into  the  program  duiing  a  USAAMRDL-sponsorcd  project,  providing  improve¬ 
ments  that  arc  directed,  principally,  toward  a  better  analytical  capability  for  studies  of 
loads,  vibrations,  and  transient  aeroclastic  behavior  of  rotor  systems. 
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DESCRIPTION  OF  MATHEMATICAL  MODEL 


RIGID  BLADE  ANALYSIS 

The  following  description  of  the  mathematical  model  applies  to  the  version  of  the 
Rotorcraft  Flight  Simulation  delivered  under  the  gust  response  contract  (Reference  2). 

The  aircraft  fuselage  is  treated  a»  a  rigid  body  that  is  subject  to  gravity  and  applied 
external  forces  and  moments  due  to: 

1.  Fuselage  aerodynamics 

2.  Main  rotor  aerodynamics 

3.  Tail  rotor  aerodynamics 

4.  Wing  aerodynamics 

5.  Elevator  aerodynamics 

6.  Fin/ruddcr  aerodynamics 

7.  Auxiliary  thrust 

8.  Weapon  recoil  force 

These  forces  and  moments,  which  act  at  input  points  of  application,  arc  transferred  to 
the  aircraft  eg  and  summed.  For  trimmed,  level,  lg  flight,  the  sum  of  the  external 
forces  and  moments  and  the  weight  muse  equal  zero.  During  maneuvers,  the  sum  is  used 
as  the  forcing  function  in  the  Euler  form  of  the  equations  of  motion  that  are  derived  in 
Reference  8. 

Aerodynamic  forces  and  moments  due  to  the  wing,  fin, 'rudder,  fuselage,  and  elevator  arc 
calculated  from  the  local  aerodynamic  environment  and  input  constants  that  provide  a 
curve  fit  to  wind  tunnei  test  data  or  to  lift,  drag,  and  moment  data  predicted  by 
analytical  methods.  The  force  and  moment  characteristics  are  generally  represented  as 
functions  of  angle  of  attack  and  Mach  number  with  corrections  applied  for  yawed  flow 
and  interference  effects.  The  forces  act  at  ccnter-of-prcssurc  locations  specified  by  the 
user. 

The  rotor  strip  theory  aerodynamic  calculations  take  into  account  an  approximate 
nonuniform  inflow  as  well  as  velocities  due  to  airspeed,  aircraft  pitch  rate,  rotor  rpm. 
arid  blade  flapping.  Only  rigid  blade  flapping  is  included  in  the  analysis.  Blade  toicc 
coefficients  are  functions  of  angle  of  attack  and  Mach  number  and  may  be  input  in 
tabular  or  curve-fit  formats.  Only  blades  with  uniform  airfoil  section  from  blade  cutout 
to  blade  tip  may  be  represented. 
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DESCRIPTION  OF  MATHEMATICAL  MODEL 


RIGID  BLADE  ANALYSIS 

The  following  description  of  the  mathematical  model  applies  to  the  version  of  the 
Rotorcraft  Flight  Simulation  delivered  under  the  gust  response  contract  (Reference  2). 

The  aircraft  fuselage  is  treated  as  a  rigid  body  that  is  subject  to  gravity  and  applied 
externa!  forces  and  moments  due  to: 

1.  Fuselage  aerodynamics 

2.  Main  rotor  aerodynamics 

3.  Tail  rotor  aerodynamics 

4.  Wing  aerodynamics 

5.  Elevator  aerodynamics 

6.  Fin/ruddcr  aerodynamics 

7.  Auxiliary  chrust 

8.  Weapon  iccoil  force 

These  forces  and  moments,  which  act  at  input  points  of  application,  arc  transferred  to 
the  aircraft  eg  and  summed.  For  trimmed,  level,  lg  flight,  the  sum  of  the  external 
forces  and  moments  r.nd  the  weight  must  equal  zero.  During  maneuvers,  the  sum  is  used 
as  the  forcing  function  in  the  Euler  form  of  the  equations  of  motion  that  arc  derived  in 
Reference  8. 

Aerodynamic  forces  and  moments  due  to  the  wing,  fin/ruddcr,  fuselage,  and  elevator  are 
calculated  from  the  local  aerodynamic  environment  and  input  constants  tiiat  provide  a 
curve  fit  to  wind  tunnel  test  data  or  to  lift,  drag,  and  moment  data  predicted  by 
analytical  methods.  The  force  and  moment  characteristics  are  generally  represented  as 
functions  of  angle  of  attack  and  Mach  number  with  corrections  applied  for  yawed  flow 
and  interference  effects.  The  forces  act  at  ccnter-of-pressurc  locations  specified  by  the 
user. 

The  rotor  strip  theory  aerodynamic  calculations  take  into  account  an  approximate 
nonuniform  inflow  as  well  as  velocities  due  to  airspeed,  aircraft  pitch  rate,  rotor  rpm. 
and  blade  flapping.  Only  rigid  blade  flapping  is  included  in  the  analysis.  Blade  force 
coefficients  arc  functions  of  angle  of  attack  and  Mach  number  and  may  be  input  in 
tabular  or  curve-fit  formats.  Only  blades  with  uniform  airfoil  section  from  blade  cutout 
to  blade  tip  may  be  represented. 


ELASTIC  ULADE  ANALYSIS 


The  analysis  developed  under  another  Army  contract  differs  from  the  rigid  blade  analysis 
primarily  in  that  treatment  of  rotor  dynamics,  hub  in-plane  motion,  shaft  windup,  and 
rotor  blade  elastic  bending  is  considered.  The  hub  motion  is  represented  bv  an  upside 
down  pendulum  with  a  torsional  spring  that  is  driven  by  the  in-plane  blade  shears  and  by 
the  pitch  and  roll  moments  as  shown  for  longitudinal  motion  in  Figure  2. 
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Figure  2.  Model  for  Longitudinal  Hub  Motion. 


The  ability  of  the  rotor  shaft  to  twist  in  response  to  the  applied  in-piane  torque  is 
simulated  by  assuming  that  at  the  bottom  of  the  shift  there  is  an  infinite  torsional 
inertia  revolving  at  a  constant  speed.  Conn-  'ting  this  large  inertia  with  the  blade  is  a 
torsionally  flexible  shaft  with  an  input  spring  rate.  The  torsional  deflection  of  the  hub 
is  then  calculated  from  the  applied  in-p!anc  blade  moment. 

In  the  rotor  analysis,  the  clastic  deformation  of  thu  rotor  blades  is  approximated  by  a 
finite  scries  of  products  of  functions  that  vary  with  radius  (mode  shapes)  and  time 
(participation  factors).  For  Z  (the  out-of-planc  defection).  Y  (the  in-plane  deflection), 
and  0  (the  torsional  deflection), 

Z(x,t)  =  Z^xJqjCt)  +  Z2(x)q2(r)  +  .  .  -  .  +  Zn(x)qn(t) 

Y(x,t)  =  Y^sJqjtt)  +  Y2(x}q2(t.  +  .  .  .  .  +  Yn(x)qn(t) 

0(x,t)  =  eaxjq^t)  +  02(x)q2(t)  +  -  -  .  .  +  0n(x)qn(t) 


where  1  ==  n 


6  and  n  is  the  number  of  modes  used. 


The  set  of  functions  [Zj(x),  Yj(.v).  ]  is  called  the  i£^’  mode  shape  of  the  rotor  blade. 

Associated  with  each  mode  shape  is  a  natural  frequency,  coj  .  Both  the  mode  shapes 
and  the  natural  frequencies  result  from  the  study  of  the  rotor  blade  free  vibration 
equations  of  motion  and  arc  used  because,  through  them,  the  solution  of  the  traced 
equations  of  motion  for  the  blade  is  vastly  simplified.  The  mode  shape  and  natural  fre¬ 
quencies  that  arc  inputs  to  the  Rotorcraft  Flight  Simulation  arc  computed  at  the  Eustis 
Directorate  by  using  the  Myklestad  program.  Other  installations  could  use  other  programs 
as  long  as  the  modes  and  frequencies  thus  calculated  are  put  into  the  Rotorcraft  Flight 
Simulation  in  accordance  with  the  User’s  Manual  instructions. 

The  displacements  and  velocities  due  to  the  hub  motion,  shaft  windup,  and  clastic  blade 
deflections  are  included  in  the  calculation  of  blade  aerodynamics,  resulting  in  a  coupled 
aeroelastic  analysis.  In  addition,  blade  loads  arc  calculated  tor  trim  and  maneuver  cases. 

A  harmonic  analysis  of  trimmed  flight  loads  is  also  provided. 


INPUT 

Input  to  the  Rotorcraft  Flight  Simulation  (Figure  3}  is  user-oriented,  i.e.,  either  physical 
measurements  taken  from  the  aircraft  or  curve  fits  to  nondimcnsional  aerodynamic  data. 
The  program  input,  generally,  includes: 

1.  Logic  control  cards 

2.  Aircraft  gross  weight;  eg  location;  fuselage  aerodynamic  center  location; 
aerodynamic  curve-fit  constants  to  fuselage  lift,  side  force,  and  drag,  and 
pitch,  roll,  and  yaw  moment  equations 

3.  Main  and  tail  rotor  physical  description,  aerodynamic  data  in  tabular  or 
curve-fit  form,  blade  frequencies  and  inode  shapes,  and  pylon  and  shaft 
stiffness  and  inertia  characteristics 

4.  Wing,  fin/rudder,  and  elevator  ccntcr-of-pressurc  areas,  aeiodynamic 
curve-fit  constants,  incidence  angles,  and  coefficients  to  approximate 
interference  effects 


5.  Auxiliary  thrust  location  (if  any),  orientation,  and  level 

6.  Control  linkages 

7.  Flight  constants  such  as  altitude,  airspeed,  and  density  and  initial  guesses 
for  trim  contro!  positions  and  aircraft  attitude 
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For  maneuvers,  die  following  additional  data  arc  required: 

8.  Weapon  location  and  orientation 

9.  Stability  and  control  augmentation  system  (SCAS)  transfer  function  constants 

10.  Description  of  the  maneuver 

11.  List  of  variables  to  be  plotted  as  functions  of  time 

The  Myklcstad  program  for  calculating  mode  shapes  and  natural  frequencies  requires  a 
complete  description  of  the  distribution  of  rotor  blade  properties,  including 

1.  Twist 

2.  Weight 

3.  Beamwise.  chordwise,  and  torsional  stiffness 

4.  Shear  center  and  eg  location 

5.  Beamwise  and  chordwise  moments  of  inertia 

Although  most  of  the  inputs  arc  easily  understood,  many  arc  nor  nearly  so  easily- 
obtained.  Experience  has  shown  that  at  least  a  week  will  be  required  to  compile  an 
operating  deck  if  data  are  readily  available.  Estimation  of  parameters  not  available  may 
require  considerable  engineering  experience,  supplementary  computation,  or  refinement  to 
obtain  correlation  with  experimental  data. 


CONTROL  LINKAGE 


PROGRAM  USE 


The  analysis  for  the  computer  program  consists  of  three  major  portions:  trim,  in  which 
the  control  positions,  aircraft  attitude,  and  blade  flapping  angles  necessary  to  achieve  a 
stabilised  flight  condition  are  calculated;  stability  analysis,  in  which  aircraft  and  rotor 
stability  derivatives  and  roots  of  the  aircraft  characteristic  equations  arc  obtained;  and 
maneuver,  in  which  the  response  of  the  aircraft  is  computed  for  nonequilibrium  flight 
conditions.  The  following  subsections  describe  uses  of  the  program  in  each  cf  these 
major  modes  of  operation. 


TRIM 

Trim  values  of  control  positions,  rotor  flapping  angles,  and  aircraft  attitude  are  computed 
for  either  level  or  climbing  flight  or  for  steady-state  turns  or  steady  "g”  pull-ups.  Aero¬ 
dynamic  forces  and  moments  from  the  rotors,  fuselage,  and  lifting  surfaces  are  compared 
with  components  of  the  aircraft  weight  (or  steady-state  . lar  accelerations  for  steady- 
state  maneuvers)  for  a  given  set  of  control  positions  and  aircraft  attitude.  Simultaneously, 
rotor  flapping  moments  and  accelerations  arc  compared  to  determine  if  they  are  consistent 
with  the  rotor  blade  root  boundary  conditions.  If  an  equilibrium  state  docs  not  result, 
the  required  changes  in  the  independent  variables  arc  computed. 

If  it  is  assumed  that  for  Z-forcc,  Z  =  Z(CP,  F/'AC,  LC,  P,  A1MR,  B1MR,  A1TR,  B1TR, 

e,*} 

where  CP  =  collective  pitch 
F/AC  =  fore-aft  cyclic 
LC  =  lateral  cyclic 
P  =  pedal 

AlMR  =  main  rotor  longitudinal  flapping 
B1MR  =  main  rotor  lateral  flapping 
A1TR  =  tail  rotor  longitudinal  flapping 
B1TR  =  tail  rotor  lateral  flapping 
0  =  fuselage  pitch  angle 
•t*  =  fuselage  roll  angle  (or  'P  =  yaw  angle) 

and  tliat  yaw  angle  is  an  input  while  roll  angle  is  a  variable,  then 


42  “  3CP  acP  +  3w5t 


AF/AC  ,  ALC  +  az  AP  +  AA1MR 


8MMR  iB1MR  +  sKfre  aA,TR  +  9tnw  AB1TR  +  1-  Ae 


+  ?§■  A*  = 


-  error  in  Z- force. 
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When  similar  equations  are  written  for  the  other  aircraft  forces  and  moments  and  the 
flapping  moments  for  both  rotors  and  arc  collected  in  matrix  form,  the  following 
equation  results: 

10  X  10  10  X  1  10  X  1 

(Partial  derivative  matrix)  (Controls  change  matrix)  =  -(Error  matrix) 

In  the  C-81  trim  procedure,  the  error  matrix  is  computed  from  the  previous  trim 
iteration;  the  partial  derivative  matrix  is  evaluated  numerically,  and  this  matrix  equation 
is  solved  to  find  changes  to  the  independent  control  variables  that  would  produce  an 
exact  trim  condition  if  the  system  were  linear.  Since  it  is  not.  more  than  one  iteration 
is  required;  but  5  to  20  iterations  arc  usually  sufficient  to  obtain  a  force  balance. 

Repetition  of  this  iterative  process  eventually  provides  values  of  the  independent  variables 
required  for  steady-state  flight.  Trim  values  of  aircraft  performance  parameters  and  rotor 
loads  are  then  printed  out.  A  multiple  case  capability  permits  incremental  changes  in 
input  parameters  for  efficient  examination  of  the  effects  of  speed,  gross  weight,  rpm,  etc. 
A  sample  of  the  trim  output  for  the  acroclastrc  version  of  the  program  is  shown  in 
Figure  4. 

This  output  suggests  the  following  uses  of  the  program: 

1.  Prediction  of  rotor  performance  and  l"ads 

2.  Prediction  of  contrc!  positions 

3.  Prediction  of  speed-power  polars 

4.  Prediction  of  aircraft  attitude  versus  speed 

5.  Definition  of  required  control  range 

6.  Definition  of  eg  limits  and  maximum  gross  weight  for  stable  flight 

7.  Prediction  of  steady-state  flight  envelope  considering  rotor  stall  and  loads, 
tail  rotor  and  stabilizing  surface  sizes,  and  power  available 

8.  Optimum  rotor  design 

9.  Analysis  of  lift  sharing  for  compound  helicopters 

10.  Prediction  of  increases  in  performance  due  to  advanced  airfoils 

11.  Evaluation  of  propulsive  force  requirements  and  effects  of  auxiliary  propulsion 

12.  Prediction  of  aircraft  turn  and  climb  performance 


13.  Design  of  tail  rotor  and  control  surface  sizes 

14.  Design  of  control  couplings  for  best  attitude  versus  speed 


STABILITY  AND  CONTROL 

A  stability  analysis  is  available  for  use  at  each  trim  condition  and  for  user-specified  times 
during  a  maneuver.  Small-perturbation  stabilility  and  control  characteristics  are  predicted 
taking  into  account  six  fully  coupled  fuselage  degrees  of  freedom  plus  pylon  and  rotor 
degrees  of  freedom.  Characteristics  for  pitch  attitude  response  to  forc-aft  cyclic,  roll 
attitude  response  to  lateral  cyclic,  and  yaw  angle  rate  response  to  pedal  input  are  also 
calculated.  A  sample  of  the  output  from  the  stability  analysis  section  of  the  program 
is  shown  in  Figure  5. 

Uses  of  the  stability  and  control  section  of  the  program  include: 

1.  Prediction  of  stability  derivatives 

2.  Orientation  and  sizing  of  stabilizing  surface 

3.  Design  of  automatic  control  system 

4.  Definition  of  aircraft  stability  limits 

5.  Prediction  of  control  power 

6.  Analysis  of  control  gearing  design 

7.  Selection  of  optimal  hub  restraint 

8.  Prediction  of  control  margin 

9.  Determination  of  effects  of  lift  sharing  on  compound  helicopter  stability 

10.  Design  of  integrated  control  system 

MANEUVER 

In  the  maneuver  section,  aircraft  response  to  nonequilibrium  flight  conditions  is  calcu¬ 
lated.  In  this  section,  the  program  is  a  true  digital  flight  simulator,  solving  a  highly 
coupled  and  nonlinear  set  of  equations  describing  the  dynamics  and  aerodynamics  of 
the  aircraft  to  predict  its  time-variant  behavior.  Maneuvers  of  the  greatest  interest  to 
most  program  users  are  movements  of  the  controls,  encounters  with  vertical  or  horizontal 
gusts,  changes  in  engine  torque  or  auxiliary  thrust,  weapon  fire,  and  activation  or 
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deactivation  of  the  automatic  stability  and  control  system.  Any  of  the  allowable  maneu¬ 
vers  may  be  used  separately  or  in  combination. 

The  user  should  be  forewarned  thar  maneuvers  can  consume  large  quantities  of  central 
processing  unit  (CPU)  time  on  a  computer.  Whereas  determination  of  a  trim  condition 
generally  requires  only  3  or  4  minutes,  simulation  of  a  3-  or  4-second  real-time  maneuver 
may  require  45  minutes  of  computer  time  on  an  IBM  360.  Also,  simulation  of,  say,  a 
pull-up  will  usually  involve  running  of  the  maneuver  several  times.  Since  flight  path  is 
not  an  input,  the  user  must  estimate  the  control  motions  required  to  produce  the  desired 
maneuver.  The  first  case  output  is  then  used  to  refine  the  input  for  the  second  case,  and 
so  foith,  until  the  desired  maneuver  is  obtained.  A  “restart”  feature  in  the  newer  ver¬ 
sions  of  the  program  allows  the  u  .cr  to  start  the  maneuver  for  the  second  case  at  any 
real-time  point  calculated  in  the  first  case,  thus  reducing  total  computer  time  (and  cost). 

A  sample  of  the  program  maneuver  output  is  shown  in  Figure  6. 

Following  are  a  few  of  the  uses  of  the  maneuver  section  of  the  Rotorcraft  Flight  Simu¬ 
lation: 

1.  Determine  aircraft  gust  response 

2.  Determine  maneuver  loads  and  biade  response 

3.  Study  minimum  distance  to  clear  an  obstacle 

4.  Study  reaction  of  the  aircraft  to  firing  on-board  weapons 

5.  Examine  large-perturbation  response  to  control  inputs 

6.  Simulate  autorotations 

7.  Determine  thrust  and  power  maneuver  limits 

8.  Simulate  tactical  maneuvers 

9.  Determine  total  aircraft  stability  by  examining  time  histories 

10.  Control  surface  sizing  for  improved  maneuverability 

11.  Study  methods  for  maximum  deceleration 

Reference  9  explains  the  use  of  the  Rotorcraft  Flight  Simulation  for  aircraft  evaluation. 
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C3NTB3U  PARTIAL  DERIVATIVE  MATRICES 

P1JN03/INCH  DR  fd.dt-rojvos/ivcm 


COLLECTIVE 

P/A  CYCLIC 

LAT  CYC.IC 

PEDAL 

<-*ORCE 

-420.1 

439,5 

-22.00 

-16.84 

j 

y-»0RCE 

-55.75 

55.70 

229,1 

3)6.0 

| 

Z-p*DRCE 

-5599. 

1171, 

-36.71 

12.69 

YA4  MOM 

191ft. 

1190. 

•247.0 

-9689. 

? 

»jt:h  mom 

1973, 

-2R51 , 

131.5 

411.4 

90„L!  mom 

-390.7 

376.5 

1548, 

1294. 

P/» 

Pj.*a  MOM 

11.02 

10.99 

7.184 

-162.1 

1 

Mc*.  L4T 

P-.A*  MOM 

-5.102 

-6,097 

-4.021 

54.12 

r.9.  E/a 

r,*’  MOM 

.0 

.0 

.0 

.0 

T.R.  L4T 

p.A*  MOM 

.0 

.0 

.0 

.0 

4.*.  P/A 

3V.3M  MOM 

13.96 

-17.43 

-1.925 

. 1 963E-01 

l 

M.9.  L4T 

PY^OV  m9M 

.9465 

.1557 

•9.053 

.6226E-02 

r.i.  p/a 

»YL0V  “DM 

.0 

.0 

.0 

.5504E-03 

t.r.  la? 

3V.jM  MOM 

.3 

.0 

.0 

-.4589E-03 

PT/SEC**2 

Q9  RAD/SEC**3  *E9  I»4 

CH 

collective 

F/A  CYCLIC 

LAT  CYCLIC 

PEDAL 

<-s0RCt 

-1.359 

1.626 

-.9159E-01 

-.624SE-01 

a 

r-roRcs 

-.2069 

.2966 

.9499 

1.172 

| 

Z-?DRCE 

-13.72 

9.344 

••  1 362 

.4706E-01 

l 

Y4*  MOM 

,1614 

.1945 

-.2192E-01 

-.7714 

\ 

»!t;m  mom 

.1123 

-.2174 

.1003E-01 

.3136E-01 

> 

40.L  ‘*'1*' 

-.1339 

.1324 

.5443 

.4552 

i 

i 

4.1.  P/A 

P.A»  MOM 

.T90PP-0? 

.7716E-02 

.5089E-02 

-.1148 

4.1.  L*T 

P.AP  MOM 

-.337QE-02 

-.4319E-02 

-.2949E-02 

.3934E-01 

f 

i.  4.  P/A 

*.*»  MOM 

.0 

.0 

.0 

.0 

: 

T.3.  I_AT 

p„Aa  mom 

.0 

,0 

•  o 

.0 

• 

4.1.  P/A 

ar-.3M  MOM 

.  1  945E+20 

•  1 152E+0? 

-135.3 

1.063 

*.4,  L*T 

pv.'jm  mom 

.31955+07 

«• 1 34UE-50 

-.9761 

.U901E-07 

T,4.  P/» 

py.OM  mom 

,167ftE-flfl 

.9371E-52 

-2.035 

-a. 391 

s 

T.4.  L*T 

ay-.D'l  MOM 

.055UE-73 

-50.25 

-135.3 

.JU71E-37 

r 

OOviV^MTIOMAL  ffTXPO  “IN 

5  M0V-0ImEV3I04Av  derivatives 

! 

collective 

P/A  CYCLIC 

•_AT  CYCLIC 

pedal 

1 

■i 

<-*04CP. 

-.5129 

.5352 

-.2885E-01 

-.2055E-01 

1 

Y-;-D3CE 

-. 59O6E-0 1 

. 6799E-0 1 

,2797 

.3956 

| 

Z-'D9CE 

-9.515 

1.030 

-.44B2E-01 

.1549E-01 

ft*  MOM 

.2129 

,1392 

-.2991E-01 

-1.017 

| 

»!T;M  mom 

,  6709 

-1.305 

.6021E-01 

•  J  883 

- 

4!)„l  M0M 

-«9a59E-0i 

, 400SE-0 1 

.1912 

.1515 

* 

4.1.  P/A 

P .A 3  MOM 

.5097E-02 

, O987E-02 

.3289E-02 

-.7420E-01 

4.1.  LAT 

p*.A3  mom 

-.7238E-03 

-.7134E-03 

-.4709E-03 

.6337E-02 

r.i.  p/a 

Pi.A3  .  MOM 

.5 

.0 

.0 

.0 

T.4.  i.AT 

P.A3  MOM 

A 

•  ^ 

.0 

.0 

.0 

*,4.  P/A 

PTuOM  MOM 

.5393E-02 

-.7990E-0C 

-,ft«19E-03 

.8989E-05 

m.4,  l*t 

py-,OM  mom 

.  1 1 08E-03 

» 1923E-00 

-.  1 06OE-O2 

.7289E-06 

I. 4.  P/A 

*Yw3M  MOu 

.0 

.0 

.2520E-06 

T.4.  LAT 

PY'.OV  MOM 

.0 

.0 

.0 

-.5373E-07 

Figure  5.  Sample  Stability  Section  Output. 
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Figure  5.  Concluded 
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VALIDITY  OF  MATHEMATICAL  MODEL 


It  must  be  emphasized  that  the  validity  of  the  correlation  presented  here,  as  well  as  the 
mathematical  model  of  any  other  aircraft  or  rotor  system,  depends  on  two  factors: 

1.  The  adequacy  of  the  techni  al  analysis  as  implemented  in  the  computer 
program. 

2.  The  validity  of  the  input  data  used  to  represent  the  physical  characteristics 
of  a  particular  aircraft. 


References  1,  2,  and  3  present  the  analysis  in  sufficient  detail  for  the  engineer  to  under¬ 
stand  the  assumptions  it  contains.  The  input  deck  also  contains  assumptions  that  arc  not 
always  clearly  visible  from  the  computed  results.  Experience  has  shown  that  in  a  large 
percentage  of  cases,  failure  of  the  Rotorcmfr  Flight  Simulation  to  obtain  a  trim  solution 
or  tc  provide  adequate  correlation  with  test  data  is  due  to  inadequate  input  data.  The 
most  frequently  made  and  most  successfully  used  suggestion  for  users  experiencing  diffi 
culty  in  running  the  program  is  “Check  your  input  data.”  In  other  v  ds.  it  takes 


+ 


Input  Data 


Math  Model 


to  represent  the  airaah,  and  the  user  must  do  his  part. 


CORRELATION  WITH  TEST  DATA 

The  correlation  data  presented  :n  this  section  demonstrate  some  of  the  virtues  and  vices 
of  the  Rotorcraft  Flight  Simulation.  They  also  reflect  the  ability  of  the  user  to  prepare 
input  for  the  program. 


Performance 


The  1/7.5  scale  CH-47C  model  rotor  correlations  presented  in  Figures  7  through  12  are 
instructive,  since  they  isolate  the  rotor  and  require  the  minimum  amount  of  input  data. 
Although  the  predicted  slopes  of  the  thrust-control  axis  angle  curves  in  Figures  7  and  10 
do  not  agree  well  with  test  results  (possibly  due  to  wake  effects),  the  correlations  in  the 
other  Figures  of  this  set  show  that  performance  prediction,  in  general,  is  good  slid  that 
prediction  of  gross  or  total  aircraft  parameters  is  likely  .o  be  acceptable.  These  curves 
also  show  the  effect  of  the  Boeing-developed  dynamic  airfoil  prediction  method  on  the 
overall  analysis. 


Figures  13  through  16  show  that  total  aircraft  level-flight  performance  prediction  is  vciv 
good  for  single-rotor  aircraft.  However,  Figures  17  and  18  show  that  for  tander.i-rotor 
helicopters,  correlation  is  good  only  for  the  forward  rotor.  The  interference  cifects  on 
the  aft  rotor  arc  not  accounted  for  in  the  analysis.  A  simple  interference  correction 
developed  by  Boeing-Vcrtol  results  in  a  much  improved  power  correlation. 


Loads 

Looking  at  the  rotor  analysis  in  greater  detail.  Figures  19  through  26  show  that  although 
pcak-to-peak  loads  arc  predicted  accurately,  the  azimuthal  variations  of  both  loads  and 
lift  coefficient  lack  agreement  with  experimental  wave  forms.  Tins  may  again  b-_  due  to 
wake  approximations  that  predict  a  more  regular  inflow  than  actually  exists. 


Maneuver 

The  predicted  maneuver  for  a  CH-47C  pull-up  in  Figure  27  agrees  well  wit.,  flight  test 
data.  In  this  figure,  the  three-per-rev  has  been  removed  from  the  test  data  by  filtering, 
while  it  has  not  been  removed  from  the  computed  results.  The  Bell  Helicopter  Company 
Model  583  pitch  response  correlation  shown  in  Figure  28  shows  a  similar  degree  of 
correlation. 


Stability  and  Control 

■■  — —d - — —  > 

The  Rotorcrr. ft  Flight  Simulation  has  been  used  extensively  with  a  high  degree  of  con¬ 
fidence  by  Bell  Helicopter  Company  to  predict  aircraft  stability  and  control  characteristics 
for  a  number  of  years.  Figures  29,  30,  and  31  show  typical  stability  and  control  corre¬ 
lation  with  flight  test  data.  The  good  agreement  between  theoretical  and  flight  test  data 
shown  in  these  figures  suggests  that  Bell’s  confidence  in  this  area  of  the  analysis  was 
justified. 


FUTURE  CORRELATION 


Several  detailed  and  systematic  correlation  efforts  using  the  Rotorcraft  Flight  Simulation 
have  been  planned  for  the  near  future.  The  first  is  the  Bell  Helicopter  Company  correla¬ 
tion  with  H  34  model  rotor  data  mentioned  in  the  following  section.  In  conjunction 
with  the  high-speed  maneuverability  flight  test  of  their  ficxicam  (hingcless)  rotor  system. 
Bell  will  also  correlate  the  C-81  predicted  aircraft  performance  and  loads  with  the  test 


A  third  correlation  will  be  conducted  by  the  Eustis  Directorate  staff  in  conjunction  with 
the  maneuverability  and  flight  loads  flight  rest  of  the  AH-1G  to  be  conducted  for  the 
Government  by  Bell.  This  effort  will  follow  the  2-ycar  flight  test  program  that  is 
scheduled  to  start  in  the  near  future. 


rHRUST  COEFFICIENT  (CT/<r) 
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ROTOR  THRUST  COEFFICIENT  (CT/«r) 


LEGEND 

O  MODEL  ROTOR  DATA 

— - BASIC  C*81 

- C-81  MftTH  RACIAL  FLOW 

-  C-81  WITH  RADIAL  FLOW 

AND  UNSTEADY  AERO 


Figuic  10.  Model  Rotor  Correlation  Witii  Elastic  Blade  Rotor 
Theory:  Rotor  Thrust  vs  Control  Axis  Angle. 
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ENGINE  SHAFT  HORSEPOWER 


ALTITUDE  =  5000  FT  M  =  ItSO  LB 
t—  C-11  RESULTS  +  SS  MS*  LOSS 
O  TEST  OATA,  USAATA  ENSINEERIW  FLtCHT  TEST 
AH-16,  PHASE  0,  PART  2,  PERFORRANCE 


Figure  15.  AH-1G  Level-Flight  Performance  Correlation. 


01  *  7511  It  (  *  11.75  SO  FT 
a  HST  tli A  - C-ll  RESULTS 


Figure  16.  AH-1J  Level-Flight  Performance  Correlation. 
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LEGEND 

O  NODEL  ROTOR  TEST  OATA 
- BASIC  C-81 

-  C-81  IITH  RADIAL  FLOI 

AND  UNSTEADY  AERO 


TEST  CONDITIONS 
CT/<r=  0.074 

M  =  0.35 
fit  -  500  FFS 

“(I)  (90)  =0‘B 


100  200 


AZIMUTH  POSITION  (+)  -  DEG 

Figure  19.  Model  Rotor  Correlation  With  Azimuthal  Cocfficient-of-Lift 
Variation  at  0.75  Radius  (G p/a  =  0.074). 


LEGEND 

O  N00EL  ROTOR  TEST  DATA 
---BASIC  C-81 

-  C-81  IITH  RADIAL  FLOI 

AND  UNSTEADY  AERO 


TEST  CONDITIONS 
C j/a  =  0.114 
H  =  0.35 
a«  =  500  FPS 

"(I)  (90)  =  °-8 


Figure  26.  UH-1D  Oscillatory  Load  Variation  With  Speed 
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Figure  27.  Correlation  of  Theory  With  CH-47  Pull-up  Maneuver. 
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Figure  28.  Bell  Model  583  Pitch  Response 
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CURRENT  DEVELOPMENT 


UNSTEADY  AERODYNAMIC  EFFECTS 

Vertcl  Division,  The  Boeing  Company  (under  Contract  DAAJ02-71-C-0045)  has  developed 
for  use  with  the  Rotorcraft  Flight  Simulation,  methodology  for  evaluating  and  predicting 
the  aerodynamic  force  and  moment  coefficients  of  an  airfoil  in  an  unsteady  flow  environ¬ 
ment.  The  method  is  especially  useful,  since  it  is  not  dependent  on  constants  that  must 
be  derived  from  wind  tunnel  tests  for  each  new  airfoil  shape.  The  analysis  has  been 
incorporated  into  the  1971  CDC  6600  version  of  the  program,  and  the  final  report  will 
be  published  in  the  near  future.  The  analysis  will  be  included  in  the  1972  version  of  the 
program  as  an  in-house  project  during  1973. 


HELICOPTER  MATHEMATICAL  MODELING  FOR  CONTROL  SYSTEMS 

Under  a  contract  to  Honeywell,  Inc.  (DAAJ02-71-C-0053),  an  improved  mathematical  tool 
tc  evaluate  future  helicopter  control  systems  is  being  developed.  The  program  will  be 
extended  to  include  the  capability  to  simulate  fluidic  stability  augmentation  systems, 
electronic  autopilot  systems,  and  control  moment  gyros.  The  contractor  will  also  provide 
a  general  optimal  control  system  design  for  demonstration  with  the  program.  The  program 
analysis  will  be  extended  to  calculate  linear  data  for  linear  analysis  and  design  and  a  trim 
procedure  that  facilitates  convergence. 


ROTORCRAFT  FLIGHT  SIMULATION  WITH  AEROELASTIC  ROTOR  REPRESENTA¬ 
TION  CORRELATION 


Bell  Helicopter  Company  (under  Contract  DAAJ02-72-C-0086)  will  establish  the  capability 
of  the  program  to  predict  helicopter  rotor  performance  and  loads  and  will  isolate  areas  of 
the  rotor  analysis  that  need  refinement.  This  will  be  accomplished  by  performing  a  corre¬ 
lation  study  in  conjunction  with  the  H-34  model  rotor  wind  tunnel  tests  to  be  conducted 
by  Sikorsky  Aircraft  under  Contract  DAAJ02-72-C-0026.  Pretest  predictions  will  be  made 
to  establish  the  ability  of  the  program  to  predict  rotor  performance  and  loads  from  design 
data.  Additional  correlation  will  be  performed  after  the  test  so  that  exact  test  points  may 
be  matched. 


IMPROVED  AERODYNAMICS  AND  AEROELASTIC  REPRESENTATION  FOR  THE 
ROTORCRAFT  FLIGHT  SIMULATION 


Under  Contract  DAAJ02-72-C-0098,  Bell  Helicopter  Company  will  improve  the  program 
in  the  areas  of  aerodynamic  representations  and  numerical  integration  techniques.  New  or 
improved  methods  for  mathematically  modeling  the  aerodynamic  forces  and  moments 
acting  on  a  helicopter  fuselage,  stabilizing  surfaces,  wings,  external  stores,  and  rotor  blades 
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will  be  developed.  Mathematical  techniques  for  solving  sets  of  differential  equations 
typical  of  those  found  in  helicopter  aeroelastic  analyses  will  also  he  studied.  A  more 
efficient  program  able  to  model  more  types  of  aircraft  in  a  larger  number  of  tactical 
maneuvers  will  result  from  this  program. 


FUTURE  WORK 


It  is  apparent  that  the  analytical  sophistication  in  at  least  two  areas  of  the  program  does 
not  match  the  refined  aeroelastic  rotor  analysis.  That  is,  the  program  is  not  of  totally 
uniform  texture.  The  first  is  in  the  representation  of  the  rotor  inflow  and  wake.  The 
simplified  equations  expressing  inflow  velocity  as  a  function  of  radius,  azimuth,  thrust, 
and  flight  path  are  not  adequate  to  predict  rotor  performance  and  loads  at  low  airspeeds. 
Rotor/rotor,  rotor/wing,  and  rotor/tail  interference  effects  all  need  to  be  predicted  more 
accurately  for  loads  and  stability  and  control. 

The  other  nonuniform  area  of  analysis  is  in  fuselage  dynamics.  The  present  rigid-body 
assumption  should  be  replaced  with  an  elastic  fuselage  representation.  Studies  of  air¬ 
craft  vibration  levels  could  be  performed  using  the  program  if  an  elastic  fuselage  repre¬ 
sentation  were  added.  This  addition  would  also  permit  more  accurate  prediction  of  tail 
rotor  angle  of  attack  that  is  necessary  for  all  tail  rotor  investigations. 
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